Formation of vortex and antivortex polarization states in ferroelectric films 
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An effort based on the kinetic model allowed by Landau-Ginzburg-Devonshire theory is presented 
here. The formation of vortex, antivortex states in ferroelectric films are reported and theoretically 
explained. Phase transition in toroidal moment is observed. Vortex-antivortex-vortex (V-AV-V) 
triplet is generated by applying inhomogeneous transverse static localized fields in x-y plane. A 
specific boundary condition is used to keep the spirit of vortex state i.e. the net polarization across 
the boundary is always zero. 
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I. INTRODUCTION 



Circularly ordered structural 
3(, Bose Einstein condensates 



phase is denoted as vortex state which is found in ferromagnets [1|, [2| , superconductors 



Kittel 



DO. 



l| has shown the possibility of circular 



4(, and in ultracold Fermi gases 
ordering due to boundary conditions of magnetic spin vectors in any ferroic material. Existence of such vortex orders 
was extended into ferroelectrics by Ginzburg et.al. 6( in 1984. 

In 2004, Naumov et al 7j defined a different order parameter toroidal moment (G) through a first principle 
calculation, where G is the vector sum of cross products between local dipole moments and corresponding position 
vectors of the dipoles. This toroidal moment passes through a phase transition during temperature variation. In lower 
temperature, polarization forms a toroidal domain configuration. They show that FE vortex states are a promisingly 
smart material that can fasten memory performance and of higher memory capacity by five orders of magnitude 
compared to memory of same size. This explanation has accelerated the ab initio explorations in the theoretical 
domain of polarization vortex structures. Gruvcrman et al 8[ in 2008 experimentally observed doughnout shape 
ferroelectric domain pattern and further numerically analyzed it by a model based on Landau- Lifshitz- Ginzburg 
equations. Existence of such vortex polar states depends on the size as well as shape of the nanodot, where strain 



and depolarizing field play crucial role[9[. Our recent observation by time dependent Landau Ginzburg theory has 



demnonstrated the formation of boundary condition induced vortex polarization states [10(. An intrinsic coupling 
condition between two transverse in plane components can reorient a head to tail 90° to a circularly symmetric vortex 
domain. 

But, magnetic, BE, superconducting systems, antivortex is also mentioned. Recently Roy et al 13( observed 



vortex-antivortex chains in permalloy rings. The difference between ferroelectric and ferromagnetic vortex dynamics, 
nearest neighbour and depolrizing field is very strong in the former system compared to the latter. But_vortex and 
antivortex pair can not form as ferroelectric toroidal moment does not break time reversal symmetry 



12|. Whereas 



in magnetic troidal moment T = l/2iV^ 
breaks both the time and space symmetry 



N 



Ti x m%, (here m, are the local magnetic dipoles at sited in N lattices) 



11 



12| has shown that 



12|. A recent ab-initio study by Prosandev et al 
external inhomogeneous field can control the toroidal moment. Here we report the results of a study of the vortex 
antivortex vortex (VAV) switching kinetics in the ultrathin ferroelectric films, which are almost two dimensional. As 
we are applying in plane electric field, the cross product does not break time as well as spatial symmetry. 

In this Letter, we explore a bilateral excersise between the two in plane polarization components under the tutelage 
of anisotropic electric fields. External electric fields arc applied in the diagonal axis of the 2D square fcrrolcctric lattice 
as shown by Roy et al 



13j a chain of the magnetization distribution around the vortex antivortex cores of the VAV 
structure in permalloy rings. The temperature is below critical tempertature. This is qualitatively different from the 
kinetics of extrinsic switching, where a conducting atomic force micriscopic tip switches ferrolectric films into a 180° 
domain. Here dipoles are switched out of plane. 

In the framework of Landau Ginzburg Devonshire phenomenology for a most symmetric case, the spatial-temporal 
evolution of the polarization components in a 2D film of the second order ferroelectric is described by these two coupled 



nonlinear parabbolic equation 



14, 



15 



dP x 
dt 

dt 



T [{AP X - BF* - CP x P 2 y ) + 6V 2 P X ] + E x 
-- r [(AP y - BP* - CPyPl) + S\7 2 P y ] + E y 



(1) 

(2) 



where A is a pheneomenological temeperature term defined for second order transition as A = Aq(T c — T). Below 
the critical temperature T c , ferroelectric phase will appear. Aq, B, C are positive constantsis. E x and E y are 
components of effective fields appearing from external electric fields and depolrizing fields acting in x and y direction 



3 

respectivcly.Hcrc T is phcnomcnological kinetic coefficient. Corresponding toroidal moment G is defined as: 

i N 

^=^E^ X ^ ( 3 ) 

i 

For a 2D film, no z component of polarization is present. Therefore, G only have z component and that will be 

= 2N Si ( X i-Pyi ~ ViPxi) 

We adopt with the following boundary conditions: 



Polarization states 


x = — L 


x = L 


y = -L 


y = L 


Vortex 


P x =0,P y =-l 


P x =0,P y =l 


P x =l,P y =0 


P X =-l,Py=0 


Anti-vortex 


P X =Q,Py=l 


P x =0,P y =-l 


P x =l,P y =Q 


P X =-l,Py=0 



TABLE I: Boundary conditions considered 



We have probed the electric field at two specific sites S\ and S2 along x = y axis in the following way: 





Si 


S 2 


E x 


-1 


1 


Ey 


1 


-1 



TABLE II: External electric field 

In order to realize phase transition in toroid moment G z , we simulate Eq.l and Eq.2 with the boundary condition 
specified in Table 1 for vortex state. We have taken total 100 sites i.e. N=100. Figure 1(a) shows net toroid moment 
G z versas a temperature dependent parameter A' = Aq(T — Tc) (A' = —A). This result is similar in spirit of the 
work by Naumov et al who has already submitted the evidence of phase transition in G z . It is evident from figure 
1(a), that at temperature T <T C (A' < 0), we observe a nonzero toroidal moment, whereas for for T > T c i.e. A' > 0, 
G z will sharply be zero. In a time domain of finite difference, we simulated Eq.l and Eq.2. We initiate A' = —2, and 
formed the vortex after 800 time steps. Then we gradually increased the value of A' by 0.1. This process continued 
until A' = 2. Figure 1(b) and 1(c) shows the equilibrium vortex and antivortex states for A' = — 1 respectively 

We are now planning to construct VAV triplet. Earlier, a combined effect of a critical depolarizing field and an 
open circuit electrical boundary condition has been mentioned, this combination can convert a out-of-plane polarized 
nanodot into a nonzero toroidal moment. Figure 2 shows a schematic external field probing diagram. Fields are in 
the plane i.e. only have two components along x and y plane. No z component is present. That is also true for 



the depolarizing field along z axis. We have taken a system of square ferroelectric lattice in a form of 2D film with 
160 points. External fields arc applied at two places SI and S2. SI and S2 are in diagonal axis (along line x=y), 
placed equal distance from center. Coordinates of SI and S2 are {15,15} and {25,25} respectively. We used the same 
boundary condition used for creating the vortex state as mentioned in table 1. Figure 3 (a), 3(b), 3(c), and 3(d) 
are showing formation of VAV triplet for time steps 1, 50, 100, and 400 respectively If the fields are withdrawn 
simultaneously, this triplet structure will converge very slowly to a single vortex. System will again capture the 
rotational symmetry which was earlier disturbed by the anisotropic electric fields. 

Fig. 4(a) and 4(b) show equilibrium P x , and P y profile in xy plane. Figure. 4(c) describes the spatial distribution 
of toroidal momcntG z in xy plane. It is clear in Fig. 4(c), there is a chain of vortex-antivortex- vortex has been formed. 
Figure 4(d) is the plot of P x and P y along the x=y diagonal line. These apparently periodic lines cut at three points 
predicting VAV locations. Both the P x , and P y attain zero polarization values at these points. It reminds us the 
superposition of two periodic waves with same frequency and amplitude but n phase difference. 

In summary, this Letter presents a method of formation of VAV triplet in thin ferroelectric films. Simple Landau 
Divonshire Ginzburg theory in a kinetic form, and special boundary conditions can produce vortex antivortex states. 
Phase transition in toroid moment has been observed. We proposed, if external fields of anisotropic nature are skillfully 
probed in small ferroelectric thin films, system, VAV triplet can be formed. Application of two electric place in close 
vicinity can produce three distinctly written zones. This can be utilized during memory application, and used as a 
source for further development of FeRAM. 
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FIG. 1: Schematic diagram shows how and where the electric field can be applied. 
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FIG. 2: (a) Phase transition in the z component of toroidal moment G where G — 
ant i vortex. 
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FIG. 3: Vortex antivortex and vortex formation at time steps (a)l, (b)50, (c)100 and (d) 400. 



